Q-ball Formation through Affleck-Dine Mechanism 
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We present the full nonlinear calculation of the formation of a Q-ball through the Affleck-Dine 
(AD) mechanism by numerical simulations. It is shown that large Q-balls are actually produced by 
the fragmentation of the condensate of a scalar field whose potential is very flat. We find that the 
typical size of a Q-ball is determined by the most developed mode of linearized fluctuations, and 
almost all the initial charges which the AD condensate carries are absorbed into the formed Q-balls, 
whose sizes and the charges depend only on the initial charge densities. 
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A Q-ball is a non-topological soliton with some con- 
served global charge in the scalar field theory 0. A Q- 
ball solution exists if the energy minimum develops at 
non-zero <fi with fixed charge Q. In terms of the effective 
potential of the field (j>, V((f>)/4> 2 takes the minimum at 
4> y£ The Q-ball naturally appears in the spectrum 
of the minimal supersymmetric standard model (MSSM) 
||. In particular, very large Q-balls could exist in the 
theory with a very flat potential Q , such as in the MSSM 
which has many flat directions, which consist of squarks 
and sleptons, carrying baryonic and/or leptonic charges 
It provides interesting attention to phenomenology 
and astrophysics 0,^,0 • 

Cosmologically, it is interesting that large Q-balls car- 
rying baryonic charge (B-ball) can be promising candi- 
dates for the dark matter of the universe, and/or the 
source for the baryogenesis P-fTpf. Moreover, they can 
explain why the energy density of baryons is as large as 
that of the dark matter (at least within a few orders of 
magnitude) . If the effective potential of the field cf) car- 
rying the baryonic charge is very flat at large <ft, as in the 
theory that the supersymmetry (SUSY) breaking occurs 
at low energy scales (gauge- mediated SUSY breaking), 
B-ball energy per unit charge decreases as the charge in- 
creases [Q. For large enough charge, such as B > 10 12 , 
the B-ball cannot decay into nucleon, and is completely 
stable, which implies that B-balls themselves can be the 
dark matter g] with charges B = 10 14 - 10 26 0, while 
baryons are created by the conventional Affleck-Dine 
(AD) mechanism. In the case of gravity-mediated SUSY 
breaking scenario, the B-balls can decay into quarks or 
nucleons, with the decay (evaporation) rate of the Q- 
ball proportional to the surface area [^2|, and if they 
decay after the electroweak phase transition, there are 
some advantages over the conventional AD baryogene- 
sis HJl(|. For example, B-balls can protect the baryon 
asymmetry from the effects of lepton violating interac- 
tions above the electroweak scale when anomalous B + L 
violation is in thermal equilibrium. Another one is that 
Q-balls with B — L charge survive the sphaleron effects 
to create the same amounts of baryon and lepton num- 



liep-ph/9909509 



bers. In either case, it is necessary for Q-balls to have 
large charges, such as Q = 10 22 - 10 28 f|,[l(J. In this 
scenario, dark matter is lightest supersymmetric particle 
(LSP) which arises from the Q-ball decay, and parame- 
ters of MSSM could be constrained by investigating the 
Q-ball cosmology p3]. Note that it is also possible to 
have stable Q-balls in the gravity-mediated SUSY break- 
ing theory depending on the details of the features of the 
hidden sector 

Those large Q-balls are expected to be created through 
AD mechanism [[l5| in the inflationary universe [p|-^0[. 
The coherent state of the AD scalar field which consists 
of some flat direction in MSSM becomes unstable and 
instabilities develop. These fluctuations grow large to 
form Q-balls. 

The formation of large Q-balls has been studied analyt- 
ically only with linear theory p|-|To|] and numerical simu- 
lations were done in one dimensional lattices Q . Both of 
them are based on the assumption that the Q-ball con- 
figuration is spherical so that we cannot really tell that 
the Q-ball configuration is actually accomplished. Re- 
cently, some aspects of the dynamics of AD scalar and 
Q-ball formation were studied in Ref. |l(| , but the whole 
evolution was not investigated, which is important for 
the investigation of the Q-ball formation. (In the con- 
text of non-relativistic Bose gas, the dynamics of drops 
of Bose-Einstein condensate, which is a non-topological 
soliton, are studied in Ref. Jl?]].) In this Rapid Commu- 
nication, we study the dynamics of a complex scalar field 
with very flat potential numerically in one, two, and three 
dimensional lattices, without assuming spherical Q-ball 
configuration. On one dimensional lattices, it is equiva- 
lent to the system independent of other two dimensions, 
so that we are observing plane-like objects. We call them 
Q-walls. Likewise, string-like objects, which we call Q- 
strings, appear on two dimensional lattices. 

First we show where the instabilities of a scalar field 
come from. To be concrete, let us assume that the com- 
plex scalar field <E> has an effective potential of the form M 



1 



V(<&) = m 4 ln 1 



m 2 



-^W + ^l 6 , (!) 



where m is the mass of the field, H is the Hubble pa- 
rameter, A is a dimensionless coupling constant, M ~ 
2.4 x 10 18 GeV is the (reduced) Planck mass, and c is 
a constant. Hereafter, we assume the matter-dominated 
universe, where H — 2/3t. This form of the potential 
arises naturally in the gauge-mediated SUSY breaking 
scenario in MSSM g. In addition to Eq.(|l|), the A-term 
(e.g., j4($ 4 + $* 4 )) is necessary for baryogenesis, since 
it makes the AD field rotate to create baryon numbers. 
Here we assume that the A-term does not crucially affect 
other terms of the potential, and take ad hoc initial con- 
ditions neglecting the A-term. The AD field thus has the 
initial charge density. Notice that the effective potential 
for the flat direction has very similar form in the gravity- 
mediated SUSY breaking scenario ||, and the features 
of the Q-ball formation are expected to be the same. 

If we note that the field is dominated by the logarith- 
mic potential and the homogeneous mode has only a real 
part (i.e., no rotational motion), neglecting the second 
and third term in Eq.(fil), the instability band is approx- 
imately estimated as Jig] 
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for large field value |$|, which is exactly the same as the 
result of Ref. || (note that there are additional insta- 
bility bands which come from the parametric resonance 
effects because of the oscillation of the homogeneous field, 
but they are subdominant effects). Therefore, the insta- 
bility band grows as time goes on, since the amplitude 
of the homogeneous mode |$| decreases as a -3 for the 
logarithmic potential. These fluctuations originate from 
the negative curvature of the logarithmic potential, which 
produces negative pressure Q . It is expected that Q-balls 
with corresponding scales are formed. 

As we mentioned, we calculate the dynamics of the 
complex scalar field on one, two, and three dimensional 
lattices. We formulate the equation of motion for the 
real and imaginary part of the field: <f> = -^{<fii + i<f>2)- 
Rescaling variables with respect to m, we have dimen- 
sionless variables 



H 



(3) 



We have exhaustedly calculated the dynamics of the 
scalar field and Q-balls for various parameters, and find 
that Q-balls are actually formed through the Affleck- 
Dine mechanism in three dimensional lattices. They 
have thick-wall profiles which are approximately spheri- 
cal, and its charge is conserved as time goes on. Here we 
show one example. Figure nl shows the configuration of 
the Q-ball at r = 10 6 on 6<P three dimensional lattices 
with A£ = f.O in the matter dominated universe. Ini- 
tial conditions are pi(0) = 2.5 x 10 6 ,<^i(0) = 0,^ 2 (0) = 




FIG. 1. Configuration of Q-balls on three dimensional 
lattice. More than 30 Q-balls are formed, and the largest one 
has the charge with Q ~ 1.96 x 10 16 



0,^(0) = 4.0 x 10 4 ,t(0) = 100. We can see more than 
30 Q-balls in the box. The charge of the largest Q-ball 
is Q ~ 1.96 X 10 16 evaluated by 



d^q 



(4) 



Though the best way to investigate the nature of Q-ball 
formation is to calculate in three dimensions, we must use 
the large-sized box to take into account low momenta ef- 
fects and the large number of lattices in order to have 
the enough resolution. Therefore, we have calculated in 
one and two dimensions, and use these results complc- 
mentarily with results in three dimensions. 

Therefore, let us now compare the evolutions of Q-balls 
in one, two, and three dimensions. As is mentioned, a Q- 
ball is a non-trivial configuration of the scalar field, which 
we can obtain by estimating the energy minimum with 
finite charge Q fixed. From Eq. (^) we can approximately 
estimate the charge of a Q-ball as 



Q = a 6 Q D ~ a?R D q ~ const 



(5) 



where q = <fii4>2 — 4>\4>2 is the charge density, and Qd is 
the charge in D dimension. Charge conservation tells us 
that Q is constant. If we assume the form of a Q-ball as 



<j>(x,t) — 0(x) exp(iut), 
the energy of a Q-ball can be calculated as 



(6) 



E= d 



V(<j>) - i^V 



J d 3 x[E grad + Vi + V 2 ] + luQ, 



(7) 
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where 



a 2 R 2 ' 



Vi ~ m log 1 



V 2 



2m 2 J 



const.. 



(8) 



When the energy takes the minimum value, the 
equipartition is achieved: E gra d ~ V\ and E gra d ~ V2. 
From these equations and the charge conservation, we 
obtain the following evolutions: 

Roca- 4 ^ D+1 \ 

woca-^WW), (9) 

which we observed have approximately the same features 
numerically [fl8|| . Therefore, the physical size R p h ys = R a 
of the Q-balls for D = 1, 2 shrinks, while remaining con- 
stant for £> = 3. Thus, stable Q-balls (three dimensional 
objects) can be formed, but Q-walls or Q-strings shrink 
because of the charge conservation. 

Figure |2| shows the power spectra of the cases of one 
dimensional lattices with the box size L = 4096 and lin- 
earized fluctuations at r = 4.5 x 10 5 and 5 x 10 5 . As 
expected, both spectra are very similar at r = 4.5 x 10 5 , 
since fluctuations have not fully developed yet. After 
they are fully developed (r = 5 x 10 5 ), the spectrum 
becomes smooth and broad because of rescattering 
(panel 2(c)). But, even at this time, the most devel- 
oped mode is the same as that of linearized fluctua- 
tions (Compare panels 2(c) and 2(d)). Therefore, we 
conclude that the size of Q-balls is determined by the 
most developed mode of linearized fluctuations when the 
amplitude of fluctuations becomes as large as that of 
the homogeneous mode, (5(f) 2 ) ~ (j> 2 . For the case of 
Fig. ^|, the typical size is k max ~ 0.04, which implies that 
Rphys ~ a(Tf)/k max ~ 7.3 x 10 3 , where t / = 5x 10 5 
is the Q-ball formation time. At this time, the ratio to 
the horizon size is ~ 10 -2 , which corresponds to the re- 
sults of Ref . H , but this value may not have important 
meaning, since the sizes of the horizon and Q-balls have 
different time evolutions. This size is consistent with 
the actual sizes appearing on three dimensional lattices, 
where the largest Q-ball has the size R p h ys ~ 1-1 x 10 4 , 
and the average size of the second to fifth largest Q-balls 
is ~ 5.2 x 10 3 . 

We also observed that almost all the initial charges car- 
ried by the condensate of AD field are absorbed into Q- 
balls formed from the fragmentations of the condensate, 
and the amplitude of the homogeneous mode is highly 
damped, which means that they carry only a small frac- 
tion of the total initial charges. In the case of Fig. [j], 
more than 95% of the charges are stored in the Q-balls. 

Actually, the charges and sizes of Q-balls depend on 
the initial value of the charge carried by the AD con- 
densate. Since the initial charge density of AD scalar is 




FIG. 2. Power spectrum of fluctuations of AD scalar field 
when the amplitude of fluctuations becomes as large as that 



of the homogeneous mode: 
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The top panels (a) 



and (c) show the full fluctuations calculated on one dimen- 
sional lattices, while the bottom panels (b) and (d) show the 
linearized fluctuations without mode mixing. 



written as q(Q) — <pi(0)ip' 2 (0), the larger the initial am- 
plitude or angular velocity of the AD condensate, the 
larger the charge stored in Q-balls. The dependences of 
the charge and the (comoving) size of the (largest) Q- 
ball on the initial charge density are shown in Fig. pi 
Here we take one dimensional lattices, so actually we ob- 
serve Q-walls. We thus estimate their conserved charges 
as Q m ax — a 3 / dxq(x) according to Eq.(||). We find 
that one large Q-ball and a few small Q-balls are formed 
in most of the cases, while one Q-ball is formed in some 
cases such that the initial charge density q(0) is relatively 
small. Open circles denote the dependence on <p' 2 (0) for 
fixed <pi(0), while solid triangles denote the dependence 
on (fii(0) for fixed y? 2 (0) with L — 256. Since both results 
show the same dependence, the only relevant variable 
which determines the charge and the size is the initial 
charge density q(0). The dashed line represents the fit- 



ted line, Q r , 



94.3(g(0)) 



1.03 



Crosses are obtained on 



L = 512 lattices. They seem to be a little larger, and 
the box size effect might be remained. Notice that those 
Q-balls with negative charge can be produced when the 
initial angular velocity of the AD condensate (or, in our 
simulations, ^(O)) l& small enough |Tgf| . 

Since the (comoving) size of the Q-balls changes as 
time goes on, as we mentioned above, the actual values 
of the size at r = 5 x 10 7 in one dimension is not so 
important. What is important is how the size depends 
on g(0). The dotted line shows the fitted line written as 
R = 9.46 x 10- 7 (<?(0)) ' 464 , though we expect the relation 
R oc ((/(O)) 1 / 3 , which means that the charge is propor- 
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FIG. 3. Dependence of charges and sizes on the initial 
charge q(0) — ipi(0)ip' 2 (0) carried by the AD condensate on 
one dimensional lattices. Open circles and solid triangles de- 
note the dependence on ifi^O) and <pi{0), respectively. Crosses 
are obtained on larger box size. R is evaluated at r = 5 x 10 r . 
Note that the relevant variable is only q(0), since both depen- 
dences lie on the same lines. 



tional to the volume as in Eq.(||). It may be one of the 
reasons for the discrepancy that the values of R for small 
charge Q-balls may have considerable error because of 
poor resolution in spatial lattices. 

In conclusion, we consider the full nonlinear equations 
of motion of the Affleck-Dine scalar field in order to 
see the formation of the Q-ball through the Affleck-Dine 
mechanism by numerical simulations. It is shown that 
large Q-balls are actually produced by the fragmentation 
of the condensate of a scalar field whose potential is very 
flat, as in the supersymmetric standard theory. 

We find that the typical size of Q-balls is determined 
by that of the most developed mode of linearized fluctu- 
ations when the amplitude of fluctuations grows as large 
as that of the homogeneous mode: (5<f> 2 ) ~ 4> 2 ■ Almost 
all the initial charges carried by the AD condensate are 
absorbed into the Q-balls formed, leaving only a small 
fraction in the form of the remaining coherently oscillat- 
ing AD condensate. We thus can constrain parameters 
of MSSM through the fraction of the baryon in Q-balls 
(/b) in the context of the Q-ball decay producing LSP 
dark matter Jl3| 

Moreover, the actual sizes and the charges stored 
within Q-balls depend on the initial charge density of the 
AD condensate, which is in good agreement with the con- 
dition of the existence of the Q-ball; that is, it exists if the 
scalar field can take non-trivial energy minimum configu- 
ration with a fixed charge. Therefore, Q-balls with huge 
charges necessary for Q-balls to be dark matter could be 



produced if the initial charge density that the AD con- 
densate carries is large enough. 

We also find that the evolution of Q-balls crucially de- 
pends on the form of their dimensions, and the stable 
Q-balls can only exist in the form of three dimensional 
objects. Smaller dimensional objects such as Q- walls and 
Q-strings shrink as the universe expands. 

Finally, we will mention the Q-axiton (the higher en- 
ergy state Q-ball) , which was studied in Ref . |]l6| . In our 
simulations, Q-balls are actually formed. The field orbits 
in the complex plane inside Q-balls are almost complete 
circles. Moreover, even if the initial AD field orbit is 
extremely oblique such that the initial angular velocity 
is very small so as to create the negatively charged Q- 
balls, circular orbits can be seen inside both positive and 
negative Q-balls |Q. It thus seems that Q-axitons may 
appear, if ever, at the very beginning of the Q-ball for- 
mation. 
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